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Transient growth of disturbances in a non-normal combustion system could trigger thermoacoustic 
instability. In this work, the effect of heat source on triggering thermoacoustic instability in a system with 
Dirichlet or Neumann boundary conditions is investigated. For this a thermoacoustic model of a premixed 
laminar flame is developed. It is formulated in state-space by expanding flow disturbances via Galerkin 
series and coupling with a flame model, thus providing a platform on which to gain insight on the system 
stability behaviors. Transient energy growth analysis is then performed by discretizing the flame model 
into distributed monopole-like sound sources. It is shown that the thermoacoustic system is non-normal 
and associated with transient energy growth. Parametric studies are then conducted to study the effects 
of (1) the flame-acoustics interaction index ,f, (2) the eigenmode number N, (3) the mean temperature 
ratio :y between pre- and after-combustion regions on the system stability behavior and non-normality. 
It is found that the maximum transient energy growth G max is increased with increased Furthermore, 
the thermoacoustic system with more eigenmodes is found to be associated with a larger G max . 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

To achieve low NO* emission and increased combustion effi¬ 
ciency, lean premixed pre-vaporized combustion technology is 
widely applied in aero-engine and gas turbines. However, they 
are more susceptible to thermoacoustic instability [1-4]. It occurs 
frequently in other combustion systems such as rocket motors, 
ramjets, boilers or furnaces. Such instability is generated due to 
the coupling between unsteady combustion process and acoustic 
disturbances present [5], As an efficient monopole-like sound 
source, unsteady combustion produces acoustic waves [6]. And 
these pressure waves [4] propagate within the combustor and par¬ 
tially reflect from boundaries to arrive back at the combustion 
zone. When unsteady heat is added in phase with the pressure 
oscillations [7], acoustical energy increases until limit cycle oscilla¬ 
tions [8] occur. Such oscillations are wanted in thermoacoustic 
engine systems [9-19], However, they are undesirable in aero¬ 
engines and gas turbines, since the oscillations may be so intensive 
that they can cause structural damage, enhanced heat transfer and 
costly mission failure to the engine systems [5,20]. 

Thermoacoustic systems have been shown to be nonlinear 
[21,22] and non-normal [23-29], The non-normality [26] is 
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characterized by non-orthogonal eigenmodes. It arises from 
unsteady heat release and/or the complex impedance boundary 
conditions. It has also been shown that in a linearly stable but 
non-normal thermoacoustic system, there can be significant 
transient energy growth of small perturbations before their even¬ 
tual decay. If the disturbance transient growth is large enough, 
thermoacoustic instability might be triggered by causing such dis¬ 
turbances to grow to amplitudes high enough to make nonlinear 
effects significant. However, different energy measures are used 
in the literature to characterize transient growth 24,30-33], 

Chu [31] argued that the disturbance energy in a viscous com¬ 
pressible flow should include the entropy fluctuations in addition 
to conventional acoustical energy, which consists of both pressure 
and velocity fluctuations. Morfey [32] suggested that the total 
disturbance energy should include the energy associated with fluc¬ 
tuations in vorticity and entropy, beside the classical acoustical 
energy. Nicoud and Poinsot [34] argued that the classical Rayleigh 
criterion is incomplete in describing the significant sources of 
fluctuating energy in a combustion system with a flow. Giauque 
et al. [35] relaxed the definition of the disturbance energy from 
Myers by considering chemical species and heat release terms. 
Initially, the energy measure used to study the transient growth 
in thermoacoustic system included the acoustical energy only 
[24,25,30]. However, there is no consensus on the definition of 
the disturbance energy. In addition to acoustical energy, entropy 
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fluctuations were introduced by Mariappan and Sujith [23] in the 
energy measure. Subramanian and Sujith [33] introduced the 
unsteady heat release fluctuations in the energy measure. However, 
the previous transient growth analysis was conducted by assuming 
there is no temperature difference between the pre- and after¬ 
combustion regions. Thus there is a need to investigate the mean 
temperature effect on quantifying the disturbance transient growth. 
Lack of such investigations partially motivated the present work. 

Experimental investigations were performed to measure the 
transient growth rate of thermoacoustic systems. Mariappan et al 
conducted such transient growth measurements on a Rijke tube 
[28 , which is a simple but widely used platform [36-39] to study 
heat-to-sound conversion. Kim and Hochgreb [41] performed sim¬ 
ilar measurements on a lean-premixed gas turbine combustor in 
Cambridge. The transient energy growth of the disturbances due 
to the non-normality cannot be predicted by the classical linear 
stability theory [21], since it provides information only about the 
long-term evolution of the eigenmodes [1,27,29], To eliminate 
the disturbance transient growth, classical linear controllers are 
not applicable [27,29,40], 

In this work, a simplified thermoacoustic system with Dirichlet 
or Neumann boundary conditions [42] is considered. A premixed 
laminar flame is confined and leads to a mean temperature jump 
from pre- to after-combustion region. In Section 2, the system 
equations are developed in state-space. The flow perturbations 
are expanded by using Galerkin series and the flame front is 
tracked by using the classical C-equation 43,44 . And unsteady 
heat release from the flame is assumed to be caused by its surface 
variation, which results from the fluctuations of the oncoming flow 
velocity. Both nonlinear and linearized flame models are discussed. 
In Section 3, transient energy growth analysis of the thermoacou¬ 
stic system is performed. The energy measure consists of both 
acoustical energy and the unsteady heat release fluctuations by 
discretizing the flame front into distributed monopole-like sound 
sources. In Section 4, the results are presented and the key findings 
are discussed. 


part. The dimensional governing equations of the thermoacoustic 
systems comprise: 


_ du dp n 
P ~dt + dx = ° 


1 tfp fdp 

?W + Pdt 


d*p 

dx 2 


(7-1) d 

c 2 dt 


(Q f mx-x f )) 


( 1 ) 

( 2 ) 


where c is the sound speed, L 0 is the length of the duct, y is the ratio 
of specific heats, S(x) is dirac delta function describing the localized 
flame or the actuator, Q f denotes the unsteady heat release from the 
flame. 

Expanding the pressure perturbation as Galerkin series 

[25,27,30] gives, 

P(M) = I^7^7r'/n(t)iA„(x) (3) 

n=l ^ nUJn 

where 


K n = 


[’ ifi(x)dx+ [ ifi(x )dx 

JO Jxr 


1/2 


(4) 


N denotes the eigenmode number. Overdot denotes the time 
derivative. The function t//„(x) are the eigen-solutions of the homo¬ 
geneous wave equation, 


1 c^p cPp 
c 2 dt 2 dx 2 


(5) 


Since the presence of the flame results in the mean temperature 
undergoing a jump from the region 0 < x < x f to x f < x ^ L 0 , substi¬ 
tuting Eq. (3) into Eq. (5) gives 


Hn- 


dx 2 


0, 0^x<x; and 



d 2 'K 

dx 2 


0, Xf<x^L 0 


( 6 ) 


Applying the pressure continuity and velocity jump across the 
flame and the Dirichlet boundary conditions at x = 0 and x = L 0 

[25,27,30], i.e. p(x, t)| x=0 = 0 and p(x, t)| x=Io = 0 leads to 


2. Description of the thermoacoustic model 

In the present work, we consider a simplified combustion sys¬ 
tem with a monopole-like flame and Dirichlet boundary conditions 
implemented as shown in Fig. 1. 

A premixed laminar conical-shaped flame acts as an acousti¬ 
cally compact heat source. It is modeled as a thin sheet and 
described by using Dirac delta function in mathematics. The flow 
variables with a tilde are the instantaneous quantities, which 
consist of a mean part denoted with an overbar and a fluctuating 
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that satisfies the condition of orthogonality, i.e. J„° ^ m (x)i// n (x) 
dx = 0 for m t 4 n. Further analysis illustrates that co„ is related 
to x f in the form of 


c 2 sin 


CO„Xf 


Cl 

-Ci sin 


cos 

co n (Lo-x f ) 

c 2 


CD n {L 0 -Xf) 


COnXf 

Cl 
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Eq. (8) can be used to characterize the frequency co n of the eigen¬ 
mode (see Fig. 2). 

Substituting Eq. (7) into Eq. (3) yields 


_ ’I'M cin ( w " x \ 

M=i Knuj n 5,111 ( c, ) 
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Replacing p(x, t) in the momentum Eq. (1) with Eq. (9) gives 


u(x,t) = 
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N UP 
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0 ^ X < Xf 
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Fig. 1. Schematic of the combustion system with an acoustically compact premixed 
flame confined and Dirichlet boundary conditions applied on both inlet and outlet. 
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Fig. 2. Variation of predicted eigenmode frequencies with temperature ratio r 2 /r, in a thermoacoustic system with Dirichlet or Neumann boundary conditions. 


Now substituting the Galerkin series expansion form of the pressure 
perturbation, i.e. Eq. (3) into Eq. (2) leads to 


-1 
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Multiplying Eq. (11) with \!/ n (x.) and integrating from x = 0 to x = L 0 
gives 
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A little algebra shows that 
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where k„ = i/i 2 (x)dxj . When Tj = T 2 (i.e. no mean tempera¬ 
ture gradient), Ci = c 2 and k„ = (1/2) 1/2 is a constant [25], However, 
when Ti ^ T 2 , k„ depends on the mode frequency co„ and the flame 
position Xf. The damping f is taken into account by assigning a 
damping parameter to each mode as 


Cn = 


1 

271 



(14) 


The instantaneous heat release rate Q/(t) from the conical-shaped 
flame [43] is given as 


Q/(t) = 0/ + Qf(t) = p„S u A f AH r = p n S u (A f +A f )AH r (15) 

where p n is the mixture density, S u is the flame speed, A f is the 
instantaneous flame surface area and AH, is the heat of reaction 
per unit mass of the mixture. The heat release rate expressed in 
the sum of a mean and perturbation can be shown as 


Of(t)_A f {t) 

Of Af 


( 16 ) 


To determine the flame surface area, the flame dynamics are cap¬ 
tured using the classical G-equation [43,44]. The oncoming flow 
to the flame is assumed to be axisymmetric and combustion is 
assumed to occur on a thin surface whose axial position at radius 
r is given by x = l(r, t) as shown in Fig. 1. The flame initiation sur¬ 
face is then described by G(x, r, t) = 0, where G(x, r, t) = x - £{r. t). 
Assuming the surface propagates itself in the direction of its normal 
n i.e. VG/jVG| into the unburnt fluid at a constant speed S u gives 
dG/dt = -u ■ VG + S U |VG| [43,44], With the assumption that the 
density change across the surface G(x, r,t) = 0 is negligible, the par¬ 
ticle velocity is given by that in the oncoming flow (u/,0) and the 
equation for ii f , the velocity distribution at the flame simplifies to 


u f - 


di(r,t) 


\ \ 


1 + 


d£(r, t) 

dr 
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If {(r, t) is known, the instantaneous flame surface area A/(t) can 
be readily calculated by 

r 2n 


A f (t) =Af+ A f (t) = + ^ (r ’ f)) ) dr (18) 


where R 0 are the radius of the duct. It can then be shown that the 
mean and fluctuating parts of the flame surface area are given as 
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( 20 ) 


where '/ = -2n\j\ - S 2 u /u 2 G . The flame front is discretized into I< 
flame elements. Each behaves like a monopole sound source. /. k is 
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the weight factor corresponding to the trapezoidal integration for¬ 
mula as given as /. k = 1/2, as l< = 1 or K, while k ¥= 1 or I(,A k = 1. 
After discretization, it can be rewritten as 


= y ' U cos (-sj k ~ ik -' 
df ^piCiK„ \ c \ J Sr 


( 21 ) 


The boundary and initial conditions are 

l{K8r,t) = 0, {(r,t=0) = {(r) (22) 


It is worth noting that the analytical method present in this work 
can also be applied to a combustion system with Neumann bound¬ 
ary conditions as discussed in Appendix A. 


3. Transient energy growth analysis of the thermoacoustic 
system 

3.1. Linearized thermoacoustic system 


The governing Eqs. (13) and (21) of the thermoacoustic system 
can be formulated in state-space form by linearizing the flame 
model via discretizing into K flame elements as 


dx _ d^ / x A \ 
dt dtVx f / 



(23) 


where the coefficient matrices L™, L TF ,L M and L BF are given in the 
Appendix B. L is a (2 N + K) x (2N + K) matrix, x is the state vector 
as given as 


x = 


" 



fClljl 

K ""!5r 

K 2 r]2 

K 2 ^’^’- 

,, Kn>1n 

"XnIJn, 

OjV 


Xa 



i,f 1 Ar,A 2 { 2 Ar-^ K Ar 

1 --- 


(24) 


The linearized thermoacoustic system can be shown to be non-nor¬ 
mal due to the fact that L (2N+K ) X(2 N+ k) does not commute with its 
adjoint [25], i.e. LL* # L*L. The system can also be shown to be non¬ 
linear if the nonlinear flame model is used. 

Note that the linearized thermoacoustic system is formulated in 
state space. This enables the energy of the flow disturbances to be 
defined clearly as discussed later and the transient growth to be 
calculated easily. In addition, the linearized governing equations 
can be used to determine the stability behaviors of the thermoa¬ 
coustic eigenmodes. This will be helpful on the design of a stable 
combustion system, of which the eigenmodes stability need to be 
predicted. 


3.2. Definition of a measure of total disturbance energy 


To characterize the transient growth of disturbances in the sys¬ 
tem, it is necessary to define a measure. The thermoacoustic sys¬ 
tem is involved with unsteady heat release Q f , which is a 
monopole-like sound source and related to the entropy fluctuation 
[31 . As discussed in the previous works [1,31,33], the unsteady 
heat release can be considered as part of the total energy contained 
in the disturbances present in the thermoacoustic system. Thus the 
energy due to the monopole-like flame can be shown as 

Ef(t) = [ L ° {Qfimx - Xf)) 2 dx = ]T(4{*<5r) 2 (25) 

Jo \ A f J k=1 

= k Sr ) 2 (26) 

k= 1 

where ^ 0 is an index denoting the flame-acoustics interaction 
intensity. It’s effect on the system stability behavior and transient 
growth is studied in the following section. 


Besides the unsteady heat release fluctuation, the system is also 
associated with acoustic waves propagating along the combustor. 
The acoustical energy [25,26,29] £ A (X/,t) per unit cross-sectional 
area consists of N eigenmodes. It can also be separated into a 
kinetic component E k , and a potential component E p . The definition 
of E tot (Xf, t) is given as 



Pi"i(x,t) j 
, P2Ul{X,t) 


dx 

dx 



!1ll\ 

K) 


= EEn(Xf,t) 


(27) 


It is apparent that £,(x/, t) is a function of both flame position X/ and 
time t. Note that the state equation of perfect gas yp = — p 2 c\ 

is used to derive Eq. (27). 

The total energy £ t ot due to the unsteady heat release and acous¬ 
tic fluctuations includes the contributions from £„ and E f . Thus £ tot 
can be expressed as 

Emit) =E a (t) + E f (t) +y f>^ r ) 2 < 28 ) 


According to Eq. (23), L is the linear stability operator of the lin¬ 
earized system. The solution to Eq. (23) is given as 


x = exp(Lf)x(0) 

The square of the norm ||x|| is given as 




; 7n 1 


(29) 


+E^^ r ] 2 =i ex P( Lf ) x (°)ii 2 ( 3 °) 


Comparing the definition of the total energy as given in Eq. (28) 
with Eq. (30) reveals that £ tot (x/,t) is related to ||x|| as 

£ t ot (x/, f) = || Qx( t) 11 2 (31) 

where Q. is a 2N x 2N diagonal matrix as given in Appendix C. 

For a given x f , the maximum total energy growth rate at time t, 
optimized over all possible initial perturbations can be expressed 
as 


G(x f , t) 


max 

Etot(X/iO) £ t ot 


>t(x/,0) 


IIQx(0)|| 2 


exp(QLQ.~ 1 f)|| 2 


(32) 


Thus the global maximum energy growth factor over a given 
time interval [0, t] is 

Gmax(Xf) = max{G(x ; , t)} (33) 

t—>(X) 

|| exp(QLQ. _1 t)|| is the 2-norm of the matrix exponential 
exp(QLQT't). To calculate G max (Xf, t), singular values decomposition 
method (SVD) is used to determine the singular values of 
explQLQf’t). Detailed information about SVD can be found in the 
previous works [26,30], 

When G max > 1, then the disturbance energy is increased for 
some time t (known as transient growth). However, if G(t) —> +oo 
as time tends to infinity, then the thermoacoustic system is 
linearly unstable. This means that the real part of the eigenvalues 
of QLQ. -1 are positive. When G max = 1, there is no transient growth 
of the disturbance energy (also known as strict dissipativity). This 
ensures that the rate of total disturbance energy change at any 
instant of time is non-positive, i.e., 


dEtot _ 1 y 

dt yp [ ln dt V°V dt \®n/ 


+ ‘ / E^[(^ r ) 2 l 

k =1 Ul 

( 34 ) 
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4. Results and discussion 


4.3. Nonlinear response of the thermoacoustic system 

The variation of the linearized flame surface area with distur¬ 
bances frequency to is shown in Fig. 3. It can be seen that the 



Fig. 3. Flame surface area variation Af(t)/Af with to and o, during one period of the 
acoustic disturbance —n 0 C ft, as the oncoming flow disturbance is given as 

u(t) = Q c (l + 0.2 sin(mt)). 


amplitude of flame surface area ratio A f (t)/Af is reduced with 
increased to and periodically varied with <j>. This means that flow 
disturbances at higher frequency pass through the flame smoothly. 
However, the flame responds strongly to the flow disturbances at 
lower frequency. 

When the nonlinear flame model is combined with the acoustic 
model, time evolution of the thermoacoustic oscillations is deter¬ 
mined, ass shown in Fig. 4(a). It can be seen that the initial distur¬ 
bances quickly grow into a nonlinear limit cycle. The frequency 
spectrum is shown in Fig. 4(b). It is apparent that there are 
harmonics present, indicating the system nonlinearity. Fig. 4(c) 
illustrates the mode-shapes at co ] and co 2 in terms of velocity 
fluctuations. Velocity antinodes are located at both ends, as 
expected. 

In order to gain insight on the generation of the limit cycle oscil¬ 
lations, Rayleigh criterion or index is generally used to characterize 
the heat-to-sound energy conversion process. For this, we define 
the phase difference dpQ between pressure fluctuation and 
unsteady heat release rate as Rayleigh index as given as 


8(co, t) = arccos 


/' 


p(coi, @)Q(ftJi, 0)d& 


\ 


p2(co u 0)d0^C a ' Q>,, 0) d&) 

(35) 


The variation of the phase 0 p q is shown in Fig. 4(d). It can be seen 
that initially the phase difference is close to 90 degree, i.e. pressure 
fluctuations and unsteady heat release are out of phase. Thus the 
initial disturbances decay (see Fig. 4(a) at t < 0.3 s). However, with 
the phase difference reduced, i.e. the pressure fluctuations and 
unsteady heat release in phase, the acoustical energy is increased 



Fig. 4. Self-sustained combustion oscillations (a) velocity limit cycle, (b) frequency spectrum, (c) velocity mode shapes at co i and co 2 , and (d) the phase difference 0 between 
unsteady heat release Q/(t) and pressure fluctuation p(x/,t) at coi with time. 
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Fig. 5. Variation of V with the flame elements number K in a thermoacoustic 
system, as x f = 0.1, L 0 = 1.0, T, = 300 K, J = 1.1 x 10 7 and N = 6. 

and finally ‘saturated’. And limit cycle thermoacoustic oscillations 
are produced. This is consistent with the well-known Rayleigh 
Criterion. 

4.2. Parametric study of transient energy growth 

Before we conduct transient energy growth analysis, it is 
important to define a measure to check the convergence with 
increased number of flame elements. The number of flame 
elements is varied in steps of 20 and the following relative change 
is calculated 

, 7/ |G max (/C) - G max (iC - 20)1 
Gmax(f^) 

Fig. 5 shows the variation of ¥' with a change in the number of 
flame elements K. It can be seen that the relative change is less 
than 0.1% for K= 350. Thus the flame elements number is chosen 
to be 400 in our following analysis. 

Parametric studies of transient energy growth are then per¬ 
formed by using the thermoacoustic model described above. There 


are 3 key parameters: (1) the flame-acoustics interaction index 
(2) the eigenmode number N, and (3) the mean temperature ratio 
y. The effects of such parameters on the transient energy growth 
are investigated by considering one parameter at a time, while oth¬ 
ers remain unchanged. 

Let’s first consider the effect of the flame-acoustics interaction 
index ./ on the system non-normality and G max . Fig. 6 shows the 
variation of the transient growth G(t) in two thermoacoustic 
systems, as the flame-acoustics interaction index ./ is set to 3 
different values. One is associated with no mean temperature 
ratio y = 1. The other is with y = 2. Note that The flame-acous¬ 
tic interaction index is randomly chosen for illustration. 

It can be seen from both systems that when the interaction 
index ./ is large, the thermoacoustic system is linearly unstable 
with the maximum transient growth G ma x —> -Poo, as denoted by 
the solid line. This means that the eigenmodes are associated with 
real positive parts. As the flame-acoustics interaction index ./ is 
decreased, the disturbances present in the system undergo tran¬ 
sient growth, i.e. G max > 1, as denoted by the dash line in the figure. 
However, further decreasing J leads to G max = 1 .0. This indicates 
the system is associated with decaying orthogonal eigenmodes. 
Similar behaviors are observed in the system with y = 2. 
However, the gradient of energy growth G(t) is different. Detailed 
study on the mean temperature effect is shown in the following 
figures. 

The eigenmode number N is then considered. Fig. 7 shows 
the variation of the transient growth G(t) in two thermoacoustic 
systems. One is involved with one eigenmode, i.e. N = 1. The 
other is with 6 eigenmodes, i.e. N = 6. It can be seen that with 
the number N of the eigenmodes increased, the maximum 
transient growth G max is slightly increased. This is most likely 
due to the increased degrees of freedom present in the 
thermoacoustic systems. The increase of G max indicates the 
intensified energy exchange between various eigenmodes. And 
the non-normality of the thermoacoustic system is increased. 
Fig. 7(b) shows the gradient of transient energy growth in the 
thermoacoustic system with 6 eigenmodes is larger than that 
with one eigenmode. This reveals that the energy conversion 
is intensified in the thermoacoustic systems with more eigen¬ 
modes present. This is most likely due to the increased degrees 
of freedom. 



Fig. 6. Time evolution of transient growth G(t) in thermoacoustic systems with (.T > 1 ) and without ( 2 T = 1) mean temperature gradient, as the flame-acoustics interaction 
index J is set to 3 different values: (a) ST = l,x/ = 0.1, Ti = 300, L 0 = 1.0 m, K = 400 and N = 6. and (b) S' = 2,x f = 0.1, T, = 300, L 0 = 1.0 m, K = 400 and N = 6. 
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(a) (b) 

Fig. 7. (a) Time evolution of the transient energy growth G(t), and (b) variation of the transient growth gradient G(t) with G(t), as x f = 0.1, Ti = 300, K = 400, L 0 = 1.0 m, and 
■/=!.! x 10 7 . 


Further analysis is conducted to gain insight on the effect of 
mean temperature ratio Two thermoacoustic systems are con¬ 
sidered. One is associated with 1 eigenmode and the other is with 
6 eigenmodes. Fig. 8(a) and (b) show the variation of G(t) with 
time, when N = 1 and N = 6 respectively. It can be seen that the 
maximum transient growth G max at = 2 is smaller than that at 
ZT — 1. This means that for a given the maximum transient 
growth G m ax is increased with decreased mean temperature ratio 
y. And the thermoacoustic system is more susceptible to trigger 
thermoacoustic instability. 

From an engineering point of view, it is important to know the 
‘safest’ heat source position, where the heat source is less suscep¬ 
tible to trigger combustion instability. To determine such posi¬ 
tions, the transient growth analysis is applied by calculating G max 
via varying x f . This is due to the fact that the maximum transient 
growth G max is shown to be nonlinearly dependent on the heat 



source position Xf/L 0 . The variation of G max with x f is shown in 
Fig. 9 at two different frequencies, it can be seen that placing the 
heat source at velocity node for example Xf = 0.5 leads to G max to 
be minimized. In other words, the ‘safest’ heat source location cor¬ 
responding to t«i is at the middle of the combustor. This finding is 
consistent with the linearized modal analysis of the thermoacou¬ 
stic system [21,38], It confirms that the model can be used to study 
the stability of the thermoacoustic system and to characterize the 
transient energy growth. 

In summary, the thermoacoustic system is shown to be non¬ 
normal and associated with disturbances transient energy growth, 
when the flame and acoustics fluctuations are interacted. And ther¬ 
moacoustic instability might be triggered by the disturbances tran¬ 
sient growth. With the flame-acoustics interaction index increased, 
the system’s stability behavior is changed from linearly stable but 
with transient growth to linearly unstable. 



(b) 


Fig. 8. Variation of transient energy growth G(t) with time in the thermoacoustic system: (a) with one eigenmode, and (b) with six eigenmodes, as x f = 0.1, as L 0 = 1 m, 
T, = 300 K, and K = 400 and J = 1.0 x 10 6 . 
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Fig. 9. Variation of maximum transient energy growth G max with heat source 
location Xf in the thermoacoustic system: (a) coi, and (b) co 2 , as L 0 = 1 m, Ti = 300 K, 
and 1< = 400 and .f = 1.0 x 10 6 . 

5. Conclusions 


dp/dx | x=0 = 0 and dp/dx\ x=Lo = 0. Similarly the basis function i p(x) 
can be shown as 


0 ^ x < Xf 

sin C0S ( M "*/A 0 Xf <x < L 0 (A1) 

Sln| c 2 ) s in(ca n (L 0 - X/ )/c 2 ) x f<X^L 0 


The eigenfrequency m„ can be determined by using the nonlin¬ 
ear relationship as 

p 2 c 2 sin(ar„X;/c 1 ) sin ( ^°^°- —— 

= p 1 c l c os(^cos( m(L °-^ 


(A.2) 


The system equation can be shown as 


d 'In | r ^k + c0 2 tl - a (y- 1 ) 
dt 2 +tn dt + n ln ~ " K„ 


where k„ is given as 
r ft 0 n '/2 


r L o 

/ 'And* 

Jo 


V2 


Qf(t) cos ( ^ 


2 




V £ 2 


(A.3) 


(A.4) 


The non-normal interaction between a premixed laminar flame 
and acoustic waves is studied in this work. The non-normality is 
characterized by transient energy growth of disturbances present 
in the thermoacoustic system. Unlike previous works, the 
unsteady heat release fluctuation is considered as part of the dis¬ 
turbance present by discretizing the flame front into distributed 
monopole-like source. The system equations are formulated in 
state-space by expanding flow disturbances via Galerkin series 
and discretizing flame front for transient energy growth analysis. 
Parametric studies are then conducted to study the effects of (1) 
the flame-acoustics interaction index J, (2) the eigenmode 
number N, and (3) the mean temperature ratio 2T between pre- 
and after-combustion regions on the system non-normality. It is 
shown that with increased eigenmode number N, the maximum 
transient growth Gm ax is slightly increased. This is most likely 
due to the increased degrees of freedom present in the system. 
In addition, the mean temperature ratio ST is found to play an 
important role in determining the maximum transient growth 
Gmax and the system stability behavior. G max is shown to decrease 
with increased 2T, as flame-acoustics interaction index ,/ is small. 
However, as ,/ is greater than a ‘critical’ value, the thermoacoustic 
system becomes linearly unstable, which is characterized by 
eigenvalues with a real positive part. In summary, the present 
work identifies 3 key parameters and studies their effects on 
determining the thermoacoustic system stability and transient 
energy growth behaviors. The reported analytical approach can 
also be applied to a thermoacoustic system with Neumann 
boundary conditions. 
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By discretizing the flame front into K elements and following the 
same analysis method introduced, the non-normality and transient 
energy growth of disturbances in the thermoacoustic system with 
Neumann boundary conditions can then be investigated. For 
simplicity, detailed analysis is neglected here. 


Appendix B. Coefficient matrices 


The coefficient matrices L™, L TF 

L ba and L bf 

are given as 
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Appendix A. Thermoacoustic system with Neumann boundary L BF 
conditions 


0 0 0 0 


The above Galerkin series expansion method can be applied to a 
thermoacoustic system with Neumann boundary conditions, i.e. 


0 0 0 0 
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Appendix C. Coefficient matrix Q 


Q is a (2N + K) x (2 N + I<) diagonal matrix as given as 
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